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Abstract 

In many physical problems it is not possible to find an exact solution. How- 
ever, when some parameter in the problem is small, one can obtain an ap- 
proximate solution by expanding in this parameter. This is the basis of per- 
turbative methods, which have been applied and developed practically in all 
areas of Physics. Unfortunately many interesting problems in Physics are of 
non-perturbative nature and it is not possible to gain insight on these problems 
only on the basis of perturbation theory: as a matter of fact it often happens 
that the perturbative series are not even convergent. 

In this paper we will describe a method which allows to obtain arbitrarily 
precise analytical approximations for the period of a classical oscillator. The 
same method is then also applied to obtain an analytical approximation to 
the spectrum of a quantum anharmonic potential by using it with the WKB 
method. In all these cases we observe exponential rates of convergence to the 
exact solutions. An application of the method to obtain a fastly convergent 
series for the Riemann zeta function is also discussed. 



1 Introduction 

In this article we review a method for the evaluation of a certain class of integrals 
which occurr in many physical problems. The method that we propose has been used 
to obtain arbitrarily precise approximations to the period of a classical oscillator, 
to the deflection angle of light by the sun and to the precession of the perihelion of 
a planet in General Relativity ^EHH], to the spectrum of a quantum potential 
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and to certain mathematical functions, such as the Riemann zeta function This 
paper is organized in three sections: in section |51 we outhne the method and explain 
its general features; in sectional we discuss different applications of the method and 
present numerical results; finally, in section 0] we draw our conclusions. 



where F = f{x±) and f{x) < F for x^ < x < x+. We also ask that v > —1 so 
that the singularities are integrable. Integrals of this kind occurr for example in the 
evaluation of the period of a classical oscillator or in the application of the WKB 
method in quantum mechanics. We wish to obtain an analytical approximation to 
lu with arbitrary precision. 

The idea behind the method that we propose is quite simple: we introduce a 
function /o(x), which depends on one or more arbitrary parameters (which we will 
call A) and define Fq = fo{x±). Although the form of /o(a;) can be chosen almost 
arbitrarily, we ask that the integral of eq. (0) with F ^ Fq and f{x) fo{x) can 
be done analytically. 

In the spirit of the Linear Delta Expansion (LDE) |S] we interpolate the original 
integral as follows: 



This equation reduces to eq. in the limit 5 = 1, however it yields a much simpler 
integral when 5 = 0. We therefore write eq. ((21) as: 



2 The method 



We consider the problem of calculating integrals of the form: 
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where we have defined 
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We can use the expansion 
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which converges uniformly for |x| < 1. 

As a result we can substitute in eq. © the series expansion of eq. © provided 
that the constraint |A(x)| < 1 is met for any x_ < x < In general, as we 
will see in the next Section, this inequality provides restrictions on the values that 
the arbitrary parameter A can take. Under these conditions the integral can be 
substituted with a family of series (each corresponding to a different A): 

~, i — n + 1 n\ 

where 

[i^o - k{x)r [A(x)]" g{x) dx . (7) 



We assume that each of the integrals defining can he evaluated analytically. 
Although we have not yet specified the form of foix), which indeed will have to 
be chosen case by case, we already know that, if all the conditions that we have 
imposed above are met we have a family of series all converging to the exact value 
of the integral I,^, after setting 6 = 1. Since the rate of convergence of the series 
will clearly depend on the parameter A, we can pick the series among all the infinite 
series representing the same integral which converges faster. In fact, although A is 
a completely arbitrary parameter, which was inserted "ad hoc" in the integral, and 
therefore the final result cannot depend upon it. When the series is truncated to 
a given finite order, we will observe a residual dependence upon A. We invoke the 
Principle of Minimal Sensitivity (PMS) [Zj to minimize, at least locally, such spurious 
dependence and thus obtain the optimal series representation of the integral: 



dX 



0, (8) 



where we have defined li^^ as the series of eq. © truncated at n = and taking 
6 = 1. We will see in the next Section that this simple procedure allows to obtain 
series representation which converge fastly. Interestingly, in general the optimal 
series obtained in this way display an exponential rate of convergence. 



3 Applications 

In this Section we consider different applications of the method described above. 
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3.1 The period of a classical oscillator 



As a first application, we now consider the problem of calculating the period of a 
unit mass moving in a potential ^(a;)^ El- The total energy IS 
conserved during the motion. The exact period of the oscillations is easily obtained 
in terms of the integral: 



T 



zdx, 



(9) 



VE-V{x) 

where x± are the inversion points, obtained by solving the equation E = V{x±). 

Clearly, the integral of eq. @ is a special case of the integral considered in the 
previous section, corresponding to choosing F = E, f{x) = V{x), g{x) = 1 and 
= —1/2. In order to test our method we consider the Duffing oscillator, which 
corresponds to the potential V{x) = 5 + t x'^- We choose the interpolating 



potential to be Vb(x) 



A(x) 



and obtain 



1 + A2 



(A-' + x^) 
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The series in eq. © converges to the exact period for A > Ao = 
since |A(x)| < 1 uniformly for such values of A and |x| < A. 

The period of the Duffing oscillator calculated to first order using © is then 



(1) 
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By setting 5 



1 and applying the PMS we obtain the optimal value of A, 
^PMS = ^^^A, which remarkably coincides with the one obtained in 8^ by using 
the LPLDE method to third order. The period corresponding to the optimal A is 

47r 



TpMS 



(12) 



v/4 + 3fiA^ ' 

and it provides an error less than 2.2% to the exact period for any value of fx and A. 
This remarkable result is sufficient to illustrate the nonperturbative nature of the 
method that we are proposing: in fact, a perturbative approach, which would rely 
on the expansion of some small natural parameter, such as fi, would only provide a 
polynomial in the parameter itself: therefore it would never be possible to reproduce 
the correct asymptotic behavior of the period in this way. 

Given that it is possible to calculate analytically all the integrals i-i/2n: we are 
able to obtain the exact series representation: 

' -l)"7r(2n- 1)!! 2A2\" ^/l ^2 ' 



n=0 
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Figure 1: Error over the period (in absolute value), defined as H 
100, for A = 10 and /i = 1 as a function of the order. The three sets are obtained by 
using the optimal value XpMS = \/3 ^A^ /2 (plus), a value A = 0.9 XpMS (triangle) 
and A = 1.1 XpMS (square). 



where 2-^1 is the hypergeometric function. Since eq. (|13|1 is essentially a power series, 
it converges exponentially to the exact result, which is precisely what we observe in 

X 100 for three different values 



LPMS- 



Figure where we plot the error H ; 
of the parameter A as a function of the order in the expansion. T^xact is the exact 
period of the Duffing oscillator which can be expressed in terms of elliptic functions. 
Corresponding to the optimal value of the parameter, XpMS = \/3 fJ.A/2, the rate 
of convergence is maximal. To the best of my knowledge eq. (|13|) corresponding 
to XpMS provides the fastest converging series representation of the period of the 
Duffing oscillator. 

We consider now the nonlinear pendulum, whose potential is given by V{9) = 
1 — cosO. By choosing the interpolating potential to be Vq{9) 



^^^9'^ we obtain 



2 cos — COS 9 
1 + A2 92-^2 



1, 



where is the amplitude of the oscillations. To first order our formula yields 

27T 6\ 2tt5 Ji(e) 



^1^(^ + 2) (i + A2)3/2 e 



(14) 
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where Ji is the Bessel function of the first kind of order 1. The optimal value of A 



in this case is given by 

and the period to first order is then 



/2Ji(9) 



TpMS = T^X 
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Despite its simphcity eq. (fT7|) provides an excellent approximation to the exact 
period over a wide range of amplitudes. 



3.2 General Relativity 

We now apply our expansion to two problems in General Relativity: the calculation 
of the deflection of the light by the Sun and the calculation of the precession of a 
planet orbiting around the Sun. We use the notation of Weinberg [Hj: 



B{r) = A-^{r) = 1 



2GM 



The angle of deflection of the light by the Sun is given by the expression 



A{r) 



r \ - B{ro) 



Br) 



-1/2 



dr 
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where ro is the closest approach. 

With the change of variable z = 1/r we obtain 



93/2 /■'/'^° 

A(p = 2 Tq' / 
Jo 



dz 



ro - r-jz2 - 2^(1 - r^z^^ 



(20) 



which is exactly in the form of eq. We introduce the potential Vo{z) = (^q+A^) 
and obtain 



A{z) 



\\z^-\)-2GM {l + rlz^) 



{rl + A2) z2 _ ^ 
^ '0 



(21) 



By performing the standard steps which are required by our method we obtain 
the optimal deflection angle to first order to be: 



A(f)pMS = -vr + 



TT 



8GM 
rovr 



(22) 



6 



Ixio" 



8xl0' 



^ 6xlo' 

-e- 
< 

4xi(r 

2xl0' 





100 200 300 400 500 

(m) 

Figure 2: Deflection angle of light obtained assuming G/c^ = 7.425 X 10-3° m/kg 
and M = 1.97 x 10^° kg as function of the closest approach rg. The solid line is 
the exact (numerical) result, the dashed line is obtained with eq. (|23j) . the dotted 
line is the post-post-Newtonian result of JHI) the dot-dashed line is the asymptotic 
result (ro oo). The vertical line marks the location of the photon sphere, where 
the deflection angle diverges. 




Exact 

-- LDE 
.... PPN2 
— Asymptotic 



corresponding to the optimal A: 

^"pMS = ~- (23) 

vr 

The surface corresponding to the closest approach for which A(f) diverges is known 
as photon sphere and for the Schwartzchild metric takes the value rg = 3GM. It is 
remarkable that eq. ((221)) despite its simplicity, is able to predict a slightly smaller 
photon sphere, corresponding to tq = SGM/vr. This feature is missed completely in 
a perturbative approach. 

In Figure 121 we compare eq. (|23)) with the exact numerical result, the post-post- 
Newtonian (PPN) result of jlUj and with the asymptotic result for very small values 
of ro (close to the photon sphere). We assume M to correspond to the physical 
mass of the Sun and G to the physical value of the gravitational constant. This 
corresponds to a strongly nonperturbative regime, where the gravitational force is 
extremely intense. The reader can judge the quality of our approximation. 

We now consider the problem of calculating the precession of the perihelion of a 
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Figure 3: Precession of the orbit of a planet assuming the values M = 1.97 x 10'^'^ kg, 
Gj(? = 7.425 X 10"^*^ m/kg and e = 0.2506 (eccentricity). The scale of reference oq 
is taken to be the semimajor axis of Mercury's orbit (ao = 5.971 x 10^° m). The 
solid line is the exact result, the dashed line is the result of eq. (|26|) and the dotted 
line is the leading term in the perturbative expansion. 



planet orbiting around the Sun. The angular precession is given by [H] 

rr+ 



AO = -27r + 2 



J r- 



y/Ajr) dr 



1 



E 1 



1 1 



(24) 



where E — , ^7 — r — -^7 — r 

are the shortest (perielia) and largest (afelia) distances from the sun. By the change 
of variable z = 1/r we can write eq. (jUJ) as 



2+ 



1 



dz 



^{z+- z){z- z_) V(l-2GM(z + z_ +Z+)) 



27r, (25) 



where z± = ^/r^. Once again the integral has the form required by our method. 
One obtains 



AO = 2-K 



Tr{3G^LM^ + a{-4:L^ + iSGLM - U7G^M^)) 

6 G M 
L 



4a(L-6GM)2Wl 



(26) 
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where a is the semimajor axis of the elhpse, given by a = (r_ + r+)/2, and L is the 
semilatus rectum of the elhpse, given by 1/L = + l/r_)/2. The optimal A is 



ApA/5 = ■ (27) 

In Figure |31 we plot the precession of the orbit calculated through the exact 
formula (solid line), through eq. (|2H|1 (dashed line) and through the leading order 
result AOo = (dotted line) ^ . Once again we find excellent agreement with 

the exact result. 



3.3 The spectrum of a quantum potential 

In 1^ the method described in this paper was applied to the calculation of the 
spectrum of an anharmonic potential within the WKB method to order h^. The 
WKB condition is 

A{E) = ^{n + 1/2) (28) 

where A(E) to order 0{fi^) is given by 

A(E) = Ji(E) J2{E) + ^ ?J3(^) . (29) 

^ ^ ^ ' A8mdE ^ ' 11520 dE^^^ ' ^ ' 

We have defined the integrals: 



X- 

XJ, 



Ji{E) = / JE-V{x)dx (30) 



^^^^^ - L ^E - vix) ' (''^ 

where x± are the classical turning points. The spectrum of the potential V{x) can 
be obtained by solving eq. H29() . The integrals appearing in the equations (|5n|) , (|^T|l 
and (|32() are of the form required by eq. (^. We can test the method with the 

2 2 4 __ 

quantum anharmonic potential V{x) = ^ — h in it was proved that the 
integrals above can be analytically approximated with very high precision with our 
method. 
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Figure 4: Error over the energy of the anharmonic oscihator with h = 1, m = 1/2, 
CO = 2 and fi = 8000 (set 1) and with h = m = uj = 1 and ^ = 4 (set 2). The boxes 
and the pluses have been obtained with our method, the triangles correspond to the 
error calculated using the analytical formula of jl2j . 



By solving eq. ()28() once that the integrals have been approximated with our 
method one obtaines an analytical formula for the spectrum of the anharmonic 
oscillator [1]: 



4 2 
1\3 / 1\3 64 



En « ei [n+-j +62 [n + -j +63 + - + • • • (33) 

where the first few coefficients are given by 

ei « 0.867146 (^) (34) 



h m 



2/3 



62 « 0.42551 I '-^-^ 1 (35) 

2 4 

63 f« -0.0466914 (36) 

64 ~ 0.030669 I ^ 1 + 0.00424238 ( ) . (37) 



In Fig. 13 we display the error over the energy defined as S = 



jp(approx) T-, (exact) 



T-,(e3:;act) 
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100 as a function of the quantum number n. The boxes have been obtained using 
our formula eq. (|33() and assuming h = 1, m = 1/2, uj = 2 and fi = 8000. In this case 
^(exact) energies of the anharmonic oscillator calculated with high precision 

in last column of Table III of The jump corresponding to n = 25 is due to the 
low precision of the last value of Table III of . The pluses and the triangles have 
been obtained using our formula eq. (pluses) and eq. (1.34) of jT^ (triangles) 
and assuming h = m = uj = 1 and = 4. In this case E^^"-'^^ are the energies of the 
anharmonic oscillator numerically calculated through a fortran code. We can easily 
appreciate that our formula provides an approximation which is several orders of 
magnitude better than the one of eq. (1-34) of We also notice that the formula 
of IE] yields a quite different asymptotic expansion in the limit of n ^ 1. We are 
not aware of expressions for the spectrum of the anharmonic oscillator similar to the 
one given by eq. H33() . 

3.4 The Riemann zeta function 

The method outlined above can be applied also to the calculation of the Riemann 
zeta function [3]. We consider the integral representation 



(-2)"-! 1 /•! log"-^ 



X 



((n) = ^ — ^ -— / — dx . (38) 

^ 2"-i - ir(n) Jo 1 + x ^ ^ 

Although eq. is not of the standard form of eq. ((3), we can write it as: 



(-2)"-_^J_ 1 log"-^ X 



C(n) = ^ — / f dx, (39) 



where A is as usual an arbitrary parameter introduced by hand. In this case A{x) 



and the condition 



x—X 



1+A 



< 1 is fullfilled provided that A > 0; one can expand 



the denominator in powers of and obtain: 

rf.l = (^)""^ ^ 1 ^ fc! (-1)^- 

' 2"-i-l^^^(l + A)^+i^^„ j!(fc-j)!^ {l+jr- ^^^^ 

Despite its appearance this series does not depend upon A, as long as A > 0. 
This means that when the sum over k is truncated to a given finite order a residual 
dependence upon A will survive: such dependence will be minimized by applying 
the PMS [7], i.e. by asking that the derivative of the partial sum with respect to A 
vanish. To lowest order one has that Ap^^^ = 2~" and the corresponding formula is 
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Figure 5: Difference |Capp(3) — C(3)| as a function of the number of terms in the sum. 



found: 



C(n) 



(2) 



n-l 



in— 1 



1 7< (1 



E 



k=0 



_2-n(A:-j) 



(41) 



We want to stress that eq. (|41|) is still an exact series representation of the 
Riemann zeta function. This simple formula yields an excellent approximation to 
the zeta function even in proximity of s = 1 where the function diverges. The rate of 
convergence of the series is greatly improved by applying the PMS to higher orders. 
In Fig. iniwe plot the difference |Capp(3) — C(3)| using eq. (|ifl]) with A = XpMS (solid 
line), A = (dashed line) and the series representation 



c(3) = E 



1 



to + 1)' 



(42) 



which corresponds to the dotted line in the plot. This last series converges quite 
slowly and a huge number of terms (of the order of 10^^) is needed to obtain the 
same accuracy that our series with XpMS reaches with just 10^ terms. We notice 
that a special case of eq. ()4U() . corresponding to A = 1, was already known in the 
literature jl.Sj . 

We can extend eq. (|in|) to the critical line, s = i + ir, and write 



c 



•1 



+ IT 



(2)- 



1 ^ (1 

J- fc=0 ^ 



1 



k 

-E 



J)!(l+i)- 



(43) 
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considered in the sum of eq. The dashed curve is obtained using A = 0.3. 



Figure 6: E = Re 



X 100, as a function of the number of terms 



In Fig. IHl we have plotted the error (in percent) over the real part of the zeta 
'£<!:Hi+!l)-C(i+-)' 



function, i.e. H = Re 



X 100, function of the number of 



terms considered in the sum of eq. H43|) . We use r = 50. The solid line corresponds 
to the formula of |13j . whereas the dashed line corresponds to using our formula, 
eq. (gni), with A = 0.3. 



In Fig. [7|we plot the difference ^^"^^ + 



eft 



+ IT 



function of A for 
The 



different values of r. In this case the series is limited to the first r terms 
optimal value of A is found close to A = 0.3. 

The two figures prove that our expansion is greatly superior to the one of P^ . 



4 Conclusions 

In this paper we have reviewed a method which allows to estimate a certain class 
of integrals with arbitrary precision. The method is based on the Linear Delta 
Expansion, i.e. on the powerful idea that a certain (unsoluble) problem can be 
interpolated with a soluble one, depending upon an arbitrary parameter and then 
performing a perturbative expansion. The principle of minimal sensitivity allows 
one to obtain results which converge quite rapidly to the correct results. It is a 
common occurrence in calculations based on Variational Perturbation Theory, like 
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the present one, that the solution to the PMS equation to high orders cannot be 
performed analyticahy. Here, however, we do not face this problem, since we have 
proved that the method converges in a whole region in the parameter space: the 
convergence of the expansion is granted as long as the parameter falls in that region. 

Although in this paper we have examined a good number of applications of this 
method to problems both in Physics and Mathematics, we feel that it can be used, 
with minor modifications, in dealing with many other problems. An extension of 
the method in this direction is currently in progress. 
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